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Abstract 

Following from recent developments in Hubalek and Kyprianou [30] 
the objective of this paper is to provide further methods for construct- 
ing new families of scale functions for spectrally negative Levy pro- 
cesses which are completely explicit. This is the result of an obser- 
vation in the aforementioned paper which permits feeding the theory 
of Bernstein functions directly into the Wiener-Hopf factorization for 
spectrally negative Levy processes. Many new, concrete examples of 
scale functions are offered although the methodology in principle de- 
livers still more explicit examples than those listed. 

1 Spectrally negative Levy processes and scale 
functions 

Let X = {Xt : t > 0} be a Levy process defined on a filtered probability 
space (fi, JF, P), where {J-'t : t > 0} is the filtration generated by X satisfying 
the usual conditions. For a; e M denote by P^; the law of X when it is started 
at X and write simply Pq = P. Accordingly we shall write E^, and E for the 
associated expectation operators. In this paper we shall assume throughout 



*University of Bath, UK 

^Centro de Investigacion en Matematicas A.C., Mexico & University of Bath, UK 



1 



that X is spectrally negative meaning here that it has no positive jumps and 
that it is not the negative of a subordinator. It is well known that the latter 
allows us to talk about the Laplace exponent ip^O) := logE[e^"^^] for 6^ > 
where in particular we have the Levy-Khintchine representation 

^Ij(9) = -ad + -a^e^ ^_ j ^^Bx _ ^ _ xei{.^^_^^)I[{dx) (1) 

2 J(-oo,0) 

where a G R, a > is the Gaussian coefficient and 11 is a measure concen- 
trated on (—00,0) satisfying ^-^(l Ax^)n((ix) < oo. The, so-called, Levy 
triple (a, a, 11) completely characterizes the process X. 

For later reference we also introduce the function $ : [0, oo) — > [0, oo) as 
the right inverse of on (0, oo) so that for all g > 

$(g) = sup{0 > : i){e) = q}. 

Note that it is straightforward to show that is a strictly convex function 
which is zero at the origin and tends to infinity at infinity and hence there 
are at most two solutions of the equation ip{9) = q. 

Suppose now we define the stopping times for each x G M 

r+ = inf{t > : Xt > x} and t~ = mi{t > : Xt < x}. 

A fiuctuation identity with a long history concerns the probability that X 
exits an interval [0, a] (where a > 0) into (a, oo) before exiting into (— oo, 0) 
when issued at x G [0, a]. In particular it is known that 



W(i){a) 



E.(e--" 1k.<.-|) = (2) 



where x G (— oo,a], q > and the function W^'^^ : M ^ [0, oo) is defined 
up to a multiplicative constant as follows. On (— oo, 0) we have W^i\x) = 
and otherwise W^''^ is a continuous function (right continuous at the origin) 
with Laplace transform 

POO 1 

/ e-'^''W^'^\x)dx = —- iov9>^q). 

Jo n^) - q 

The functions {W^^") : g > 0} are known as scale functions and for conve- 
nience and consistency we write W in place of W^^^ . Identity ([2|) exemplifies 
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the relation between scale functions for g = and the classical ruin problem. 
Indeed setting g = we have that "^xiTo < r^) = 1 — W{x)/W{a) and so 
assuming that ip'{0+) > and taking limits as a t oo it is possible to deduce 
that ly(oo)"^ = 'ip'{0+) and hence 

F.,iT^ < oo) = l-ij'{0+)W{x). 

It is in this context of ruin theory that they make their earliest appearance 
in the works of [64], [H] and then later either explicitly or implicitly in the 
work of [24] and [36l [37], [60] and [53j, the real value of scale functions as a 
class with which one may express a whole range of fluctuation identities for 
spectrally negative Levy processes became apparent in the work of [HI [T5] 
and [Zl [8] and an ensemble of subsequent articles; see for example [43j, [3j, 
[iTliSlgglEO], [40], [20] and [IHEIIES]. Moreover with the advent of 
these new fluctuation identities and a better understanding of the analytical 
properties of the function W^'^^ came the possibility of revisiting and solving 
a number of classical and modern problems from applied probability, but now 
with the underlying source of randomness being a general spectrally negative 
Levy processes. For example, in the theory of optimal stopping [21 [3] and 
[39) . in the theory of optimal control [4], [HI [51] and [46], in the theory of 
queuing and storage models [2^ and [5], in the theory of branching processes 
[TT] . [31] and [44], in the theory of insurance risk and ruin |17J, ^35], [34] , 
[20] . [40] . [51] and [46], in the theory of credit risk [29] and [42] and in the 
theory of fragmentation [38] . 

Although scale functions are now firmly embedded within the theory of 
spectrally negative Levy processes and their applications, and although there 
is a reasonable understanding of how they behave analytically (see for exam- 
ple the summary in [39]), one of their main failings from a practical point 
of view until recently is that there are a limited number of concrete exam- 
ples. None the less, Hubalek and Kyprianou [30] have been able to describe 
a parametric family of scale functions in explicit detail using a mathemati- 
cal trick in which a spectrally negative Levy process is constructed having a 
particular pre-determined Wiener-Hopf factorization. 

In this paper, we make use of the aforementioned trick and combine it 
with classical and recent developments in the potential analysis of subordi- 
nators. Stemming from the latter, the theory of special Bernstein functions 
feed directly into the theory of scale functions and brings about the exis- 
tence of pairs of spectrally negative Levy processes whose scale functions are 
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conjugate to one another in an appropriate sense. 

Following an analytical exposition of the latter concepts, we are able to 
offer completely new explicit examples of conjugate pairs of scale functions. 

According to Lemma 8.4 in ^39j, for any g > 0, the g-scale function of a 
spectrally negative Levy process can be obtained via the 0-scale function of 
a spectrally negative Levy process that drifts to oo, whose law is obtained 
by an exponential change of measure of the law of the original Levy process. 
From the theoretical point of view this allows us to restrict our study to 
only the case g = 0. Nevertheless, from the point of view of applications 
this represents a restriction given that to make such construction tractable 
one needs to know explicitly the bivariate Laplace exponent associated to 
the descending ladder process, its potential measure and the inverse of the 
Laplace exponent ip, which in general is an open problem. 

In the sequel we will assume, unless otherwise stated, that the spectrally 
negative Levy process does not drift to —oo. The reason for this is that the 
0-scale function of a Levy processes drifting to — oo can be deduced from the 
0-scale function corresponding to the associated Levy process conditioned to 
drift to oo. This will be made precise in Section [4] below. 

2 Descending ladder height and parent processes 

The principal idea for generating new examples of scale functions in Hubalek 
and Kyprianou [30], which we borrow here, relies on constructing a spectrally 
negative Levy process around a given possibly killed subordinator which plays 
the role of the descending ladder height process. For the convenience of the 
reader we devote a little time in this section reminding the reader of the 
meaning of a descending ladder height process and give the result of ^30j in 
detail. 

It is straightforward to show that the process X — X. '■= {^t — 2^t • ^ ^ 0}, 
where := mis<tXs, is a strong Markov process with state space [0, oo). 
Following standard theory of Markov local times (cf. Chapter IV of [7]), it 
is possible to construct a local time at zero for X — X which we henceforth 
refer to as L = {Lt : t > 0}. Its inverse process, L^^ := {L^^ '■ t > 0} where 
L^^ = inf{s > : > t}, is a (possibly killed) subordinator. Sampling 
X at we recover the points of minima of X. If we define Hf = X^-i 
when L^^ < oo, with Hf = oo otherwise, then it is known that the process 
H = {Ht : t > 0} is a (possibly killed) subordinator. The latter is known as 
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the descending ladder height process. Moreover, if T is the jump measure of 
H then 



/•oo 

T(x, oo) = e*(°)" / e-*(°)"n(-cx), -u)du for x > 0, 

J X 

see for example [63] • Further, the subordinator has a drift component if and 
only if cr > in which case the drift is necessarily equal to cr^/2. The killing 
rate of if is given by the constant E(Xi)VO = ^'(0+)V0. Observe that in the 
particular case where $(0) = the jump measure of H has a non-increasing 
density. 

The starting point for the relationship between the descending ladder 
height process and scale functions is given by the Wiener-Hopf factorization. 
In 'Laplace form', for spectrally negative Levy processes, this can be written 

as 

m = io-m)m, o>o, (3) 

where (j){6) = — logE(e"^^^). In the special case that $(0) = 0, that is to 
say, the process X does not drift to —oo or equivalently that ^''(0+) > 0, it 
can be shown that the scale function W describes the potential measure of 
H. In other words 

/"OO 

/ dt ■ F{Ht e dx) = W{dx) for x > (4) 
Jo 

or equivalently 



oo -I 

e-^'^Widx) = —— for e>0. (5) 



We henceforth restrict the discussion to the case that $(0) = 0, see however 
Section m 

The following theorem, taken from Hubalek and Kyprianou [30], shows 
how one may identify a spectrally negative Levy process X (called the parent 
process) for a given descending ladder height process H. An equivalent result 
has been obtained in Proposition 7 in [10] . Note that the version of the 
theorem we present here constructs the parent process such that it does not 
drift to — oo. However, this is not a necessary restriction in the original 
formulation of this result. 
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Theorem 1. Suppose that H is a subordinator, killed at rate k, > 0, with 
jump measure which is absolutely continuous with non-increasing density, say 
V, and drift d. Then there exists a spectrally negative Levy process X that 
does not drift to —oo, henceforth referred to as the parent process, whose 
descending ladder height process is the process H . The Levy triple (a, a, 11) of 
the parent process is uniquely identified as follows. The Gaussian coefficient 
is given by 

a = VM. 

The Levy measure is given by 

n(— oo, —x) = v{x), for X > 0. 

Finally 

a = / xll{dx) — K. 

J(-oo,-l) 

The Laplace exponent of the parent process is also given by 

m = om 

for 9 >0 where (piO) = — logE(e^^^^). The parent process oscillates or drifts 
to oo according to whether 0(0) = or > 0. 

Note that when describing parent processes later on in this text, for prac- 
tical reasons we shall prefer to specify the triple (a, 11, ip) instead of (a, a, 11). 
However both triples provide an equivalent amount of information. 



3 Special and conjugate scale functions 

In this section we introduce the notion of a special Bernstein functions and 
special subordinators and use the latter to justify the existence of pairs of so 
called conjugate scale functions which have a particular analytical structure. 
We refer the reader to the lecture notes of Song and Vondracek [58] and the 
books of Berg and Forst [6] and Jacob [3^ for a more complete account of 
the theory of Bernstein functions and their application in potential analysis. 

Recall that the class of Bernstein functions coincides precisely with the 
class of Laplace exponents of possibly killed subordinators. That is to say, a 
general Bernstein function takes the form 

(P(e) = n + d9+[ (1 - e-^^)T(c/a;) for ^ > (6) 
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where k > 0, d > and T is a measure concentrated on (0, oo) such that 



Definition 1. Suppose that 0(6*) is a Bernstein function, then it is called a 

special Bernstein function if 

9 > 0, (7) 



where (f)*{0) is another Bernstein function. Accordingly a possibly killed 
subordinator is called a special subordinator if its Laplace exponent is a 
special Bernstein function. 

Note that if is a special Bernstein function with representation as given 
in ((71) then one says that the Bernstein function 0* is conjugate to 0. More- 
over it is apparent from its definition that 0* is a special Bernstein function 
and is conjugate to 0*. In [26] and |59j it is shown that a sufficient con- 
dition for to be a special subordinator is that T(x, cxo) is log-convex on 
(0,oo). 

For conjugate pairs of special Bernstein functions and 0* we shall write 
in addition to ((7l) 

^*(e) = K* + d*e+ [ il-e-^'')T*{dx), 9>0, (8) 

where necessarily T* is a measure concentrated on (0, oo) satisfying J^^ oo)(-'-^ 
x)T*{dx) < oo. One may express the triple (/t*,d*,T*) in terms of related 
quantities coming from the conjugate 0. Indeed it is known that 




K > 

-1 



(0,oo) 



xT(dx)j K = 



(9) 



and 

d > or T(0,oo) = oo 

(k + T(0,oo))"^ d = and T(0,oo) < oo. 

Which implies in particular that hk* = = dd*. In order to describe the 
measure T* let us denote by W{dx) the potential measure of 0. (This choice 
of notation will of course prove to be no coincidence) . Then we have that W 
necessarily satisfies 

W{dx) = d*Soidx) + {k* + T*(x, oo)}dx for x > 0. 
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Naturally, if W* is the potential measure of 0* then we may equally describe 
it in terms of (fi:,d, T). A proof of these facts and other interesting results 
can be found in [57j. Besides, it can be easily shown that a necessary and 
sufficient condition for a function to be a special Bernstein function is that 
its potential measure has a density on (0, oo) which is non-increasing and 
integrable in the neighborhood of the origin; see e.g. [9j Corollaries 1 and 2 
for a proof of this fact. 

We are interested in constructing a parent process whose descending lad- 
der height process is a special subordinator. The first parts of the following 
theorem and corollary are now evident given the discussion in the current 
and previous sections. 

Theorem 2. For conjugate special Bernstein functions and 0* satisfying 
^ and ^ respectively where T is absolutely continuous with non-increasing 
density, there exists a spectrally negative Levy process that does not drift to 
—oo, whose Laplace exponent is described by 

m = ^ = ^'^W fore>0 (10) 

and whose scale function is a concave function and is given by 

W{x) = d* + K*x+ [ T*{y,oo)dy. (11) 
Jo 

Conversely, if ip is the Laplace exponent of a spectrally negative Levy process 
that does not drift to — oo and its associated scale function, W, is a concave 
function then there exists a pair of conjugate Bernstein functions and 0*, 
satisfying ^ and ^ respectively, where T is absolutely continuous with 
non-increasing density, and such that the relations [TW and [Tl\) hold. 



Proof. Only the second part of the theorem needs a proof. Let ip and W be 
as described in the second part of the statement of the theorem and let 
denote the Laplace exponent of the descending ladder height subordinator 
associated to the spectrally negative Levy process with Laplace exponent ip. 
The latter and former functions are related via ^ with $(0) = 0. By an 
integration by parts it follows that 

e / e-'^''W(x)dx = W(0+) + / W'(y)e-^ydy = — -, ^ > 0, 
Jo Jo n^) 
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where W denotes the first derivative of W, which exists almost everywhere 
because of the concavity of W. This implies in particular that in (0, oo) 
the potential measure of the descending ladder height subordinator has a 
density which is non-increasing as well as an atom at zero of size Vr(0+). 
Furthermore, again by an integration by parts it follows that 

POO 

/ e-^''W{x)dx 
Jo 

POO 

= eW{0) + e / W'{y)e~^ydy 
Jo 

POO 

= W'{oo) + eW{0)+ / {l-e-'^y)d{-W'{y)), 

Jo 

So, that the function 0* defined by 

is a special Bernstein function conjugated to (j). □ 

Note that the proof of the converse statement in Theorem [2] says a little 
more than is claimed. Indeed we have that the potential measure of the 
subordinator with Laplace exponent (/>*, in (0, oo), admits a density which 
is decreasing and convex. To see this, one should recall that the potential 
measure associated to 0* has a decreasing density which is given by the tail 
Levy measure of 0, and since iIj{9) = 9(f)(9), 6* > it follows by an integration 
by parts that the Levy measure of has a decreasing density. 

The assumptions of the previous theorem require only that the Levy and 
potential measures associated to have a non-increasing density in (0, oo), 
respectively; this condition on the potential measure is equivalent to the 
existence of 0*. If in addition it is assumed that the potential density be a 
convex function, in light of the representation (fTTi ). we can interchange the 
roles of and 0*, respectively, in the previous Theorem. The key issue to 
this additional assumption is that it ensures the absolute continuity of T* 
with a non-increasing density. We thus have the following Corollary. 

Corollary 1. // conjugate special Bernstein functions and 0* exist sat- 
isfying ^ and ^ such that both T and T* are absolutely continuous with 
non-increasing densities, then there exist a pair of scale functions W and 



(12) 

^ > 0. 
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W* , such that W is concave, its first derivative is a convex function, [TH) is 
satisfi,ed, and 

W*{x)=d + Kx+ T{y,oo)dy (13) 
Jo 

whose parent processes are given by ^^) and 



^*W = :i77JT = ^<^*W- (14) 



Conversely, if ip is the Laplace exponent of a spectrally negative Levy process 
such that its associated scale function W is a concave function whose first 
derivative is a convex function, then there exists a pair of conjugate Bernstein 
functions and 0* and a function such that ip* is the Laplace exponent 
of a spectrally negative Levy process, and ip (respectively i^i*) is related to 
(p* (respectively to (p) by equation flU) (respectively by [I4\)) and the scale 
functions associated to ip and ip* satisfy equations fTl\) and ( f73|j . respectively. 

Proof. The first part of the proof is a simple consequence of Theorem [21 The 
converse part follows from the calculations used in the proof of Theorem [21 
Indeed, if ip and W are as described in the converse statement of the corollary 
then it follows from Theorem [2] that there exists a pair of conjugate Bernstein 
functions (p and (p* and by hypotheses and equation ( [T2ll that the respective 
Levy measure of (p and (p* have non-increasing densities. That the function 
Ip* defined by ip*{9) = 9(p*{9), > 0, is the Laplace exponent of a spectrally 
negative Levy process is then a consequence of Theorem [H It follows from 
Theorem [2] and the uniqueness of Laplace transform that the scale functions 
associated to ip and ip* have the claimed properties. □ 

There are a number of remarks that are worth making regarding the above 
theorem in the setting of conjugate special Bernstein functions. 

1. For obvious reasons we shall henceforth refer to the scale functions 
identified in ([6| and ([8j) as special scale functions. 

2. Similarly, when W and W* exist then we refer to them as conjugate 
(special) scale functions and their respective parent processes are called 
conjugate parent processes. This conjugation can be seen by noting that 
thanks to ([7|). 

W*W*{dx) = dx. 
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3. It is known (cf. Chapter 8 of [39]) that any scale function has a dis- 
continuity at the origin if and only if the parent process has paths of 
bounded variation. That is to say, in the representation ([T]), a = and 

0) k|n((ix) < oo. Taking account of the description in Theorem [T] 
one may easily deduce that in terms of the descending ladder height 
process this is equivalent with the fact that T(0, oo) < oo and d = 
which, within the context of Theorem [2l is equivalent to the case that 
d* > in dni) as predicted by the general theory. 

4. Another known general property of scale functions is that 1^(0+) < oo 
if and only if, in the representation ^,a>Oora = andIl{—oo, 0) < 
oo. See e.g. Exercise 8.5 in [39]. The latter conditions in terms of the 
descending ladder height process are respectively equivalent to d > 
or d = and dT {0+) / dx < oo. 



4 Tilting and parent processes drifting to — oo 

In this section we present two methods for which, given a scale function and 
associated parent process, it is possible to construct further examples of scale 
functions by appealing to two procedures. 

The first method relies on the following facts concerning translating the 
argument of a given Bernstein function. 

Lemma 1. Let (f) be a special Bernstein function with representation given 
by Then for any (3 > the function 4>f3{0) = (f){6 + f3), 6 > 0, is also a 
special Bernstein function with killing term Kp = (f){l3), drift term d^ = d and 
Levy measure T p{dx) = e~^^T{dx)^ x > 0. Its associated potential measure, 
Wp, has a decreasing density in (0, oo) such that Wp{dx) = e~'^^W'{x)dx, 
a; > 0, where W denotes the density of the potential measure associated to <p. 
Moreover, let (jf and (j)*^, denote the conjugate Bernstein functions of (f) and 
(f)f3, respectively. Then the following identity 

/"OO 

(l)*^(0) = (l)*(0 + /3)-(l)*(l3)+l3 [l-e-^'')e-^''W'{x)dx, ^ > 0, (15) 

Jo 

holds. 

Proof. We will use the theory of change of mesure for Levy processes; for 
background on this topic see e.g. [55] Section 33. Assume that under P, 
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H is di subordinator with Laplace exponent 0. It follows from the indepen- 
dence and homogeneity of the increments of H that the process Mf = e~^^^ , 
t > 0, is a submartingale in the natural filtration generated by H. Moreover, 
limj^o+E (e~^-^*) = 1. It follows that the function hi3{x) := e~^^,x G [0, oo], 
under the assumption that e~°° = 0, is an excessive function for the semi- 
group of the process H. We denote by P^^^ the law of the /i-transform of P via 
the excessive function hp, that is P'^'^^ is the unique measure over the space 
of right continuous left-limited with lifetime paths such that 

p(/3) = MiP over s < t), for t > 0. 



It is easily verified that under P^^^ the law of H is that of a subordinator killed 
at rate Indeed, as P^^-' is locally absolutely continuous with respect to 

P it follows that under P*^'^^ H has non-decreasing paths and for s, t > 

for every 6 > Q and A G a{Hu,u < t). Which proves at once that under 
p(/3)^ if has independent and homogeneous increments, its Laplace exponent 
is given by = 4>{P + ■), and its killing term is thus 0/3(0) = 

Moreover, the equality 



poo 

= <P{(3) + ed+ / (l-e^^^)e-^^T(cix), ^ > 0, 



justifies the description of the characteristics of claimed in the Lemma [TJ 
Furthermore, the potential of the process H under P*^^^ is given by 



(/'oo \ poo poo 

e~f^^n{H,edx}dt^ = d*6o{dx) + e-'^'-'W {x)dx, 



for X > 0, owing to Fubini's Theorem. Given that the function x i— >■ 
e~^^W\x), X > is a decreasing function, it follows that 0/? is a special 
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Bernstein function. Furthermore, following the description in Section [3] the 
characteristics of its conjugate, (p*^ are given as follows: its killing term is 
n*p = 0, as k/^ = (t){(3) > 0, the tail of its Levy measure is given by 

T}{x,oo) = e-^''W'{x), x>0, 

and its drift equals d* as 

(j)*0{0) 1 1 
d^i = lim — - — = lim — — — = hm '* 



/3 e^oo 9 e^oo0(/? + ^) e- 

in the obvius notation. The description in f fTSll follows by bare-hands calcu- 
lations using the latter facts. □ 

Note in particular that if T has a non-increasing density then so does T/j. 
Moreover, if W convex (equivalently T* has a non-increasing density) then 
W'p is convex (equivalently has a non-increasing density). These facts 
lead us to the following Lemma. 

Lemma 2. // conjugate special Bernstein functions (j) and (p* exist satisfying 
^ and such that both T and T* are absolutely continuous with non- 
increasing densities, then there exist conjugate parent processes with Laplace 
exponents 

Md) = ^Md) and r^io) = o<PW. > 0. 

whose respective scale functions are given by 

Wi3{x) =d* + [ e'^yT*{y, oo)dy, x > 0. 
Jo 



and 



/X / POO \ 
I^J e-^'T{dz)j dy. (16) 



using obvious notation. 



The second method builds on the latter to construct examples of scale 
functions whose parent processes drift to —oo. 

Suppose that is a Bernstein function such that 0(0) = 0, its associated 
Levy measure has a decreasing density and let /? > 0. Theorem [T], as stated 
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in its more general form in [30], says that there exists a parent process, say 
X, that drifts to — oo such that its Laplace exponent ip can be factorized as 

It follows that is a convex function and ?/'(0) = = iIj{I3), so that {3 is 
the largest positive solution to the equation ip{6) = 0. Now, let Wjs be the 
0-scale function of the spectrally negative Levy process, say Xp, with Laplace 
exponent ipf3{0) := ip{6 + for 6* > 0. It is known that the Levy process Xp 
is obtained by an exponential change of measure and can be seen as the Levy 
process X conditioned to drift to oo, see chapter VII in [7]. Thus the Laplace 
exponent ipp can be factorized as ilJpiO) = OcppiO), for 6* > 0, where, as before, 
0/3(-) := + ■). It follows from Lemma 8.4 in [39], that the 0-scale function 
of the process with Laplace exponent ip is related to Wp by 

W{x) = e^'^Wpix), x>0. 

The above considerations thus lead to the following result which allows for 
the construction of a second parent process and associated scale function over 
and above the pair described in Theorem [21 

Lemma 3. Suppose that (p is a special Bernstein function satisfying ^ such 
that T is absolutely continuous with non-increasing density and k = 0. Fix 
/9 > 0. Then there exists a parent process with Laplace exponent 

whose associated scale function is given by 

W{x) = d*e^" + e^" / e-^yT*{y, oo)dy, x > 0, 
Jo 

where we have used our usual notation. 

Now, when we assume furthermore that the potential density associated 
to is a decreasing and convex function, or equivalently that the Levy mea- 
sure of (p* has a decreasing density, there are three choices for a conjugate 
parent process. The first, is the one appearing in Lemma [2] with Laplace 
exponent given by ipp{9) = 6(f)*p{6), for 6* > 0, and its scale function, 
is described in equation (fT6l) . This parent process drifts to oo. The scale 
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functions W/s and Wp are conjugated in the sense described in Remark 2 in 
Section [3l which implies that 

die-'^^'Wix)) * dW^{x) = dx, x> 0. 

The second and third candidate parent process, based on 4>*{9) are the one 
drifting to —oo, constructed using the formulation above, and the one which 
drifts to oo described in Corollary [T], accordingly as 0*(O) = or 0*(O) > 0, 
respectively. For these parent processes the respective scale functions are 
such that 

e-P=^dx - 1^ (^~^''^(^)) * ^ (e-^^lV*(x)) , if (j)*{0) =0, 

^ ^~ \d{e-^'=W{x))*{e-'^''W*{dx)), if0*(O)>O, ^~ 

Remark 1. Observe that the construction explained in this section can be 
performed as soon as there exists a (3 > such that the function 6 i— (f){(3 + 6) 
is a special Bernstein function with a non-increasing Levy density. Which in 
view of the calculations carried in the proof of Lemma [T] could occur without 
(p being a special Bernstein function in itself. 



5 Complete scale functions 

We begin by introducing the notion of a complete Bernstein function with a 
view to constructing scale functions whose parent processes are derived from 
descending ladder height processes with Laplace exponents which belong to 
the class of complete Bernstein functions. 

Definition 2. A function is called complete Bernstein function if there 
exists an auxiliary Bernstein function rj such that 

0(^) =e^ f e-%(x)rfa;. (17) 

i(0,oo) 

It is well known that a complete Bernstein function is necessarily a special 
Bernstein function (cf. [3^) and in addition, its conjugate is also a complete 
Bernstein function. Moreover, from the same reference one finds that a ne- 
cessary and sufficient condition for to be complete Bernstein is that T 
satisfies for a; > 

T{dx) = \[ e~^y-f{dy)\ dx 

U{0,oo) J 
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where J^^ -^^ ll{dy) + J^^^-j ^lidy) < oo. Equivalently T has a completely 
monotone density. Another necessary and sufficient condition is that the 
potential measure associated to has a density on (0, oo) which is completely 
monotone, this is a result due to Kingman |33] and Hawkes [27j. The class of 
infinitely divisible laws and subordinators related to this type of Bernstein 
functions has been extensively studied by several authors, see e.g. [I2], [62] , 
[54] . [T8] . [45] and the references therein. 

Since necessarily T is absolutely continuous with a completely monotone 
density, it follows that any subordinator whose Laplace exponent is a com- 
plete Bernstein function may be used in conjunction with Corollary [H The 
following result is now a straightforward application of the latter and the fact 
that from (fTTl) . any Bernstein function rj has a Laplace transform {9"^ / (f){9))~^ 
where is complete Bernstein. 

Corollary 2. Let i] be any Bernstein function and suppose that (f> is the 
complete Bernstein function associated with the latter via the relation ( flTlJ . 
Write (f)* for the conjugate of (j) and rj* for the Bernstein function associated 
with (j)* via the related Then 

W{x) = i]*{x)andW*{x) = ri{x), x>0. 

are conjugate scale functions with conjugate parent processes whose Laplace 
exponents are given by 

m = ^ = Om and riO) = ^ = e> O. 

We conclude this section with some remarks about the above corollary. 

1. For notational consistency we call the pair W and W* (conjugate) com- 
plete scale functions and their respective parent processes (conjugate) 
complete parent processes 

2. In essence only part of the above corollary is of practical value. That is 
to say, any given Bernstein function rj is a scale function whose parent 
process is the spectrally negative Levy process whose Laplace exponent 
is given by ip*{0) = 9^/(f){9) where is given by ffTTl) . 

3. Observe that any given completely monotone function say p : [0, oo[—y 
[0, oo], such that p{s)ds < oo, may be seen as the potential density 
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of a subordinator. This is due to the fact that any such function is 
log-convex, owing to Holder's inequality, and then by a result due to 
Hirsch [28] there exists a subordinator, say H, whose potential measure 
admits p as a density in (0, cxd), see also ^59j for a recent proof of the 
latter fact. Let be the Laplace exponent of H, this is such that the 
tail of its Levy measure is a completely monotone function that we will 
denote by p*. Thus, p*{s)ds < oo and 



d + - + y^ e-'yp*{y)dy, > 0. 



It follows from Corollary [2] that there exists a spectrally negative Levy 
process with Laplace exponent ip{9) = 9(l){9) for 6* > 0, its associated 
scale function is a Bernstein function, and can be represented as 

W{x) = d* + n*x+ ra-e-'^y)^^, a;>0; 

Jo y 

under the assumption that p admits the representation 

p(x) = K* + e~''y-f{dy), x>0, 



where k* > and 7 is a measure over (0, 00) such that 

I '-/{dy) + / < 00, equivalently / p{t)dt < 00. 

Jo Ji y Jo 

For the respective conjugates we have that 

POO n 

c^*{e) = K* + d*e+ / ^^{dy), 9>o, 
Jo Q + y 

POO 

(T* = y2d*, a* = j ye-y-i{dy)+p{l) - k\ 
Jo 

POO 

n*(-oo, -x)= / ?/e"^^7((iy), x>0, 
Jo 



and therefore, 



PX 

W*{x) = d + KX+ / p*iy)dy, x>0. 
Jo 
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Finally, owing to Remark 2 in Section [3l if d = = d*, the functions p 
and p* are related by the Volterra-type equation 

ry ry 

/ p{x)p*{y — x)dx = 1 = / p*{x)p{y — x)dx, y > 0- 
Jo Jo 

4. Another source of examples comes from the observation that if is a 
complete Bernstein function then the function 0(^) = {ip{l/9))~^ , for 
> 0, also is a Bernstein function. This assertion is easily proved using 
that the potential density associated to is a completely monotone 
function. Thus, given a pair of complete conjugate Bernstein functions 
ip, (f* the functions 0, as defined above, and 0* constructed analogously, 
form also a conjugate pair of complete Bernstein functions. So that 
having knowledge of the parent processes and scale functions associated 
to (p and (f*, respectively, one can construct a new family of conjugate 
parent processes and scale functions. 



6 Concrete examples 

The previous sections have essentially consisted of re-dressing the theory 
of Bernstein functions in the language of scale functions. In this section 
we justify the value of the previous exposition by offering a large cache of 
remarkably explicit examples. 

Example 1 

Our first example will not be as interesting as other examples and has been 
included principally for the purpose of illustrating how the theory works in 
the context of an 'old favorite'. 

Consider the, apparently trivial, Bernstein function 

(j)(e) = K + de 

where d, k > 0. This is the Laplace exponent of the ladder height process 
with a parent process consisting of a Brownian motion with coefficient \/2d 
drift at rate k, 

^p(e) = Ke + d^^ ^ > 0. 
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This process is a diffusion and gives us an example where its scale function as 
a spectrally negative Levy process coincides precisely with its scale function 
as a diffusion. We have by a simple Laplace inversion 

This tells us that d* = 0, W'{x) = + T*{x,oo) = d'^e-'"^/'^ and hence 
K* = 0. Regarding the conjugate parent process, it is clear that 

n*(-cx),-x) = 4e"""/^ x>0. 
d"' 

and so the latter has a compound Poisson jump structure with negative 
exponentially distributed jumps having parameter n/d and arrival rate n/d"^. 
Since k* = d* = the conjugate parent process is an oscillating process with 
no Gaussian component and hence one may write down from this information 
directly 

After a little algebra one finds that this coincides with the expected expres- 
sion given by ip*{9) = 9'^/4>{9). It follows from Theorem [2] that the scale 
function associated to ip* is given by 

W*{x) = d + Kx, x>0. 

So that W* is an ultimately linear scale function, thus its associated potential 
density is ultimately constant and so it is closely related to the potential 
measures appearing in Section 3 in [59]. Finally by the Continuity Theorem 
for Laplace transforms it follows that the the scale functions are continuous 
in the Levy triple (k, d, T). For this reason as k the scale functions 
W and W* converge towards the scale functions x t— > x/d, and x d, for 
X > 0, respectively. Which provide two more examples of ultimately linear 
scale functions. 



Example 2 

Let /3, c> 0, z/ > and A G (0, 1). We claim that 

cf3eT{u + (3e) 



^ > 0, 
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is a Bernstein function where T{u) denotes the usual Gamma function with 
parameter m > 0. In order to determine the triple (k, d, T) in ([6]) associated 
with (p let us recall that the Beta function is related to the Gamma function 
by the following formula, for a,b > 



B{a,b) := / x^~\l-xy-^dx 
Jo 



r(a)r(6) 

r(a + b) ■ 



We thus have that 

cf3e 



B{p9 + iy,X), 9>0. 



Then making a change of variable in the expression for the Beta function we 
reach the identity 

POO 

/ e-^^e-^^/^(l-e-^/^)'"'d^. (18) 
Jo 



6 r(A) 

This means that k = d = and 

T(x, oo) = c^^^ (e-/^ - 1) " , X > 0. 

It is clear from the above expression that T has a density which is monotone 
decreasing. 

In order to determine the potential measure associated to this subordi- 
nator observe the following elementary identity: 

e _T{u + [3e + x)T{i-x) u + (3e ^^^^ 



4)i9) cpT{u + pe + i) r(i-A)' 
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Therefore, we have that 



c/?r(i - 


A) 


A 




c(3T{l - 


A) 


A 




c/3r(i - 


A) 


A 





X 



oo 



(1 - ^"'') . . ' ^A+W ^ (20) 



e 

c(3r{l - A) 7o ^ ' (e- -1) 

^ r(^/ + A) A /-^ . _ e"(^-")/^ 

c/3r(z/) + c/?2r(l - A) io ^ ' ^ (e-//^ - 1)'+^ 

This shows that is a special Bernstein function whose conjugate 0* has 
triplet {K*,d*,T*) (cf. (IHD) given by n* = T{u + \) /cpT{u), d* = and 



cP^r{l - A) (e-//5 - 1)' 



T*((ia;) = — a+T'^-^ 



Note that T* has a decreasing density. Referring back to Theorem [2] and 
Corollary [T] we may now say the following. 

There exists an oscillating spectrally negative Levy process with Laplace 
exponent 

which has no Gaussian component and its Levy measure, 11, satisfies 

n(-oo,-x) = 

c(A - 1) e-'l^+WiJ 

a r(Ar' ' 
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and the associated scale function is given by 



Y(u + X) ( poo ^ (,z{l^L')/f3 



n*(— oo, — x) = — — x+T for X > 



for X > 0. 

There exists a spectrally negative Levy process which drifts to oo or 
oscillates according to whether > or z/ = 0, with Laplace exponent 

which has no Gaussian component and whose Levy measure, 11*, satisfies 

c/?2r(l - A) (e-//3 - 1)^ 
and the associated scale function is given by 

rx z(v+X-l)/p 

W\x) = ^^^^ [e^'^-lf-'dz. 

Note in the special case that z/ = and /3 = c = 1 we have the two 
conjugate scale functions 

and 

Another special case worthy of remark is the case where = 1 = /3, 
c = r(l + A). The Laplace exponent -0* takes the form 

where a = 1 + A G (1, 2). It was shown in Chaumont et al. [I6j that this is 
the Laplace exponent of a spectrally negative Levy process. The associated 
scale function, also identified in the latter paper, may simply be written 

W*{x) = {l-e-''Y~\ x>0. 
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The previous calculations can be used to provide an example of the technique 
developed in Section [H For c, z/ > 0, A g]0, 1[, let be the function defined 
by means of 

It follows from the previous discussion, that the function $ defined by 

is the Laplace exponent of a subordinator such that its Levy measure has 
a decreasing density and its potential measure has a decreasing and convex 
density. According to our discussion in Section [4] and the previous facts 
it follows that "^{0) is the Laplace exponent of a spectrally negative Levy 
process that drifts to — oo, and its scale function is given by 

„^ux rx 

W{x) = — — / e-'^l - e-'f-^dz, a; > 0. 

r(A) Jo 

This is due to the fact that the Levy process with Laplace exponent \l/ con- 
ditioned to drift to oo has Laplace exponent 

Moreover, the conjugate parent and ladder height process have Laplace ex- 
ponent given by 

^ ' v{v + \ + 6y ^ ' T{v + \ + ey - ' 

so that their characteristics were discussed before, and the corresponding 
scale function is described in equation (|2T1 ). taking [3=1. The Laplace 
exponent \l/ can be taught as the one of the parent process with Laplace 
exponent ip* , as in equation f l22l) . conditioned to drift to —oo. The particular 
case where z/ = 1 = /3, has been studied in [TB]. 

Example 3 

Let 0<a</3<l, a, 6>0 and be the Bernstein function defined by 

(j)[e) = a^^-° + 6^^, e>o. 
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That is, in the case where a < /3 < 1, is the Laplace exponent of a subordi- 
nator which is obtained as the sum of two independent stable subordinators 
one of parameter (3 — a and the other of parameter /3, respectively, so that 
the killing and drift term of are both equal to 0, and its Levy measure is 
given by 

If a = /5 < 1 then is the Laplace exponent of a stable subordinator 
killed at rate a; when a < /3 = 1, is the Laplace exponent of a stable 
subordinator with positive drift 6; and finally in the case where a = 1 = P, (f) 
is simply the Laplace exponent of a pure drift subordinator killed at rate a. 
The latter case will be excluded because it has been discussed in Example 1. 
In all cases the underlying Levy measure has a density which is completely 
monotone, and thus its potential density, or equivalently the density of the 
associated scale function W, is completely monotone. 

In the remainder of this example and subsequent examples we shall make 
heavy use of the two parameter Mittag-Leffler function defined by 



where a,P > 0. The latter function can be identified via a pseudo-Laplace 
transform. Namely, for A G M and 3^(6') > A^/" — 7, 



POO 

/ e-^"e-^"x'^-^E„,;3(Aa;")rfa; 
Jo 



(0 + 7) 



a-fS 



^ + 7)° - A 

The associated scale function to can now be identified via 

W\x) = ^x^^^E^,^ (-ax°/6) , X > 0, (23) 

which is a completely monotone function because it is the product of the 
completely monotone functions x^~^ and Eq,,/3(— x"), and the later is com- 
pletely monotone because it is the composition of the completely monotone 
function t i— Ea^p^—t) for t >0, see [56j, with the Bernstein function x". So, 
the function 

^{6) = ^0(0) = a^^-"+i + be^+\ e>o, 
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is the Laplace exponent of a spectrally negative Levy process. We shall elab- 
orate on the features of the aforementioned parent process below according 
to three parameter regimes. 

In the case a < P < 1, the parent process oscillates and is obtained by 
adding two independent spectrally negative stable processes with stability 
index (3 + 1 and 1 + (3 — a, respectively. The scale function associated to it 
is given by 

W{x) = \ ! t'^'^E^,f3{-ar/b)dt, x>0. 
Jo 

The associated conjugates are given by 
and 

The subordinator with Laplace exponent 0* has zero killing and drift terms 
and its Levy measure is obtained by taking the derivative of the expression 
in (i23ll . By Theorem [T] the spectrally negative Levy process with Laplace ex- 
ponent ip*, oscillates, has unbounded variation, has zero linear and Gaussian 
terms, and its Levy measure is obtained by derivating twice the expression 
in (I23l). 

In the case, a = /3 < 1, the Laplace exponent takes the form 

^(6) = e(f){e) = ae + be^+\ 

The latter is the Laplace exponent of a oscillating spectrally negative a-stable 
process with stability index a = (1 + /?), and positive drift with rate a. The 
scale function can be implicitly found in Furrer [25] and it takes the form 

W{x) = -(1 - E^,i(-ax^/6)), X > 0. 
The respective conjugates are given by 



and 



W*{x) = ax+ ^^^_^^ x'''', x>0. 
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The conjugate subordinator and spectrally negative Levy process, can be 
described using similar reasoning to that of the previous parameter regime 
and we omit the details. One may mention here that by letting a | the 
Continuity Theorem for Laplace transforms tells us that for the case 0(6') = 
b9^, the associated ip is the Laplace exponent of a spectrally negative stable 
process with stability parameter 1 + /3, and its scale function is given by 

^(^' = W(rw/' " - °- 

The associated conjugates are given by 

(f)*(e) = b-^e^-^, ^*{9) = r^^'-^, ^ > o, 

and 

So that 0*, respectively tp* , corresponds to a stable subordinator of parameter 
1 — /3, zero killing and drift terms; respectively, to a oscillating spectrally 
negative stable Levy process with stability index 2-/3, and so its Levy 
measure is given by 

Lastly, in the case «</? = !, the Laplace exponent of the parent process, 
tp, is associated to the addition of a spectrally negative stable process with 
stability index 2 — a plus an independent continuous Levy process with no 
drift and Gaussian coefficient equal to \/b. Its associated scale function is 
given by 

Wix) = 7 / E„i(-at"/6)dt, X > 0. 
Jo 

The respective conjugates are given by 



and 
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The characteristics of conjugate parent process can be determined as in the 
other two parameter regimes but what is different is that it is a process with 
bounded variation since H^*(0) = d = 6 > 0. 

To complete this example, observe that the change of measure introduced 
in Lemma [T] allows us to deal with the Bernstein function 

(f)(e) = k{e + mf-"" + b{e + m)^ , e>o, 

where m > is a fixed parameter. In this case we get that there exists a 
spectrally negative Levy process whose Laplace exponent is given by 

ilj{e) = ke{e + m)^-" + be{e + m)'^, e>o, 

and its associated scale function is given by 

1 r 

W{x) = - e-"'Hf^-^Eo,,p{-ar/b)dt, x>0. 
" Jo 

The respective conjugates can be obtained explicitly but we omit the details 
given that the expressions found are too involved. 

As in the Example 2, the degree of generality on which this example 
has been developed allows us to provide another example of the technique 
developed in Section [H For, m,a,b > 0, < a < P < 1, there exists a parent 
process drifting to —oo and with Laplace exponent 

^{e) = {e-m) {ae^-'' + be^^) , ^ > o. 

It follows from the previous calculations that the scale function associated 
to the parent process with Laplace exponent \E' is given by 

pmx rx 

W{x) = ^ / e-"'H'^-^Ea,i3{-ar/b)dt, x>0. 
b Jo 

Finally observe that the function defined in (l23l) is a completely monotone 
function, which would have allowed us to present this example performing 
the construction indicated in the Remark 3 in Section [5l 

Example 4 

This example builds on the work of Boxma and Cohen [13] and the gener- 
alization thereof by Abate and Whitt [I]. In the latter, a parent process is 
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considered whose Laplace exponent satisfies 

^{e) = e e>o. 

This corresponds to a process which has a linear unit drift minus a compound 
Poisson process of rate A with jumps whose distribution, F, has Laplace 
transform 1 — 6'/(/i + V0){1 + VO). Let ri{x) = e^eTk{y/x). The tail of the 
jump distribution satisfies 



F{x,oo) = yY^J ^^^^"^ - f^vixf^ )) 
which in the case fi = 1 should be interpreted in the limiting sense so that 
F{x,oo) = {2x+ l)r]{x) -2J-, x>0. 

V TT 

In both cases, the distribution also has mean and hence the mean of the 
Levy process is 1 — A/yU which is assumed to be strictly positive (so that the 
process drifts to infinity). We may thus identify the characteristics of the 
Laplace exponent, 0, of the descending ladder height process. Specifically 
K = 1/fi, d = and T{dx) = XF{x, oo)dx. The scale function associated 
with is given in [T] by 

where 

Conveniently it is shown in [1] that F{x, 00) is completely monotone which 
makes a complete Bernstein function. Moreover we automatically we get 
the existence of a conjugate complete scale function 

Jo [Jy 1-/^ J 

(with the obvious interpretation when /i = 1). The associated parent process 
has Laplace exponent 




rio) 



ifi- \)+9 + {l+ fx)y/9' 



28 



It is not difficult to show that VF(0+) = 1 and hence d* = 1 showing that 
the conjugate parent process has a Gaussian component with coefficient -\/2. 
As usual, n*(— cxD, —x) may be computed by considering W'{x) — W'{+oo) 
which happens, in this case, to be a rather cumbersome expression but, none 
the less, explicit. 



Example 5 

Let a e (0, 1), K > 0, d > 0, c > and (p be the Bernstein function defined 

by 

0(e) ^K + d9 + cr , 9>0. (24) 

(Although the case d = has been treated in Example 3 we include it here 
again because a different approach is proposed.) That is, (f) is the Laplace 
exponent of a subordinator which is a a-stable subordinator with drift d 
killed at rate k. In this case 



V'(e) ^Ke + de^ + ce^+'^, e > o, 



so that the parent Levy process associated to <j) drifts to cg and is the sum 
of an independent Gaussian process with drift and a (1 + Q;)-stable process. 
Given that the Levy measure of has a density which is completely mono- 
tone it follows that its potential measure in (0, oo) has a density which is 
completely monotone, which we will next describe. Let X be an a-stable 
subordinator, be an exponential random variable of parameter and as- 
sume that X and are independent. Observe that the random variable 
Z — eJ'^Xc -\- de^ has a density and its Laplace transform has the form 

It follows that the scale function, W , associated to has density given by 
W'{x) = hz{x)/K where hz is the density of Z and can be written using 

that Xc c^/'^Xi in terms of the density of Xi, say Pa, or in terms of the 
Mittag-Leffler function, following if d > or d = 0, respectively, as follows 

(-1) ^''^7^^ ^ ^^x"-^E.,. (-fx"), ifd = 0, 
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owing to the fact that the Laplace transform of Xf°, is given in terms of the 
Mittag-Lefiler function. An expression for the density Pa in series form can 
be found in equation (14.31) of Sato [55]. It follows that the scale function 
associated to the spectrally negative Levy process with Laplace exponent ip 
is given by P(Z < x)/k. That is to say 

for X > 0. The conjugate Laplace exponent 0* has zero drift and killing 
terms and its Levy measure is described by T*(x, oo) = W'{x), x > 0, the 
conjugate parent process oscillates and has Laplace exponent given by 

^ ^ K + d9 + c9'^ 
Finally, the conjugate scale function is given by 



c 



W*{x) = d + Kx + — -x'-°, x>0. 

r(2 — a) 

In the sequel fix 7, c > 0, d, k > 0, a G (0, 1). We would like to determine 
the scale function associated to the Bernstein function 

(j)^(e) = K + de + c{e + -fy -c-f"", e>o. (26) 

To that end assume first that k > d7 + 07°^, and then observe that 

0,(^)= 0(^ + 7), 9>0, 



where has the form (12411 with killing term k — d'j — 07". Note that 0^ is 
the Laplace exponent of the subordinator with killing term k, drift term d 
and Levy measure 



cae 



--yx 



T^(dx) = — —r—dx, X > 0. 

^ r(l-a)a;i+" 

Note the density given above is the product of completely monotone functions 
and hence is itself completely monotone thus making 0^ a complete Bernstein 
function. The associated scale function, W^, has a density in (0, cx)) which 
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is given by W!^{x) = e '^^W'{x). It follows from (l25ll and the fact that 
d; = limeioo l/<Pjid) = 0, k* = hmg^o = lime^o = 0. 

i;/;^""'"'!'°"'E.,,(-^!/")<i!/, ifd = 0, 

for x > 0. These also gives some insight about the form of the scale function 
for any value of 7,c > 0. Indeed, integrating directly in the case d > 0, 
and term by term in the case d = 0, it is easily seen that for any k, c > 0, 
and 7 > the function 

^>.(?#)pfe, ifd>0, 




has Laplace transform 1/0^(6*), with (p^ as defined in (p6ll . So, for any a G 
(0, 1), K, c, 7 > 0, d > 0, the scale function associated to 0^ is given by (l27ll . 

On account of the continuity theorem of Laplace transforms, the scale 
function W is continuous in the Levy triple (k, d, T). For this reason, when 
d = 0, taking further k | 0, we recover from (l27l) 



W^{x) = - [ e-^^y"-iE„,,(7V)% 
c Jo 



which is one of the scale functions found in [30] . 

The conjugate scale function W* is easily computed taking account of the 
Levy triple associated with 0^ to be 



Observe that taking c = A/a in the definition of 0^, in (l26l) . and making 
a tend to 0, we get that 

lim ( K + d^ + -((0 + 7)"-7") ) = A; + d^ + Alog(l + - 1 , ^ > 0. 
\ a / V 7/ 

So, by the continuity of Laplace transforms it follows that the scale function 
converges as a — to the scale function corresponding to the Bernstein 
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function k + d6 + X log (^1 + - j , which is studied in the case k = = d, 

7 = 1, in the forthcoming Example 7. 

A different approach to this example would start with the Bernstein func- 
tion in (l25l) . which by the second remark following Corollary [2] we know is 
the scale function associated to some spectrally negative Levy process, and 
then determine its associated parent process, ladder height process and con- 
jugates. 



Example 6 

Another example belonging to a related family of Bernstein functions to those 
of the previous example starts again with calculations found in [30]. We take 

0(^) = . + A(i-(^y), e>o, 

where k, A > 0. For this (j) we also have d = and 

T{dx) = ^x'^-^e-^^c/x, a; > 0, 

where z/ G (0, 1) and 7 > 0. Note that the assumption on u ensures that 
T has a non-increasing density and hence (j) may be used as a descending 
ladder height process. This Bernstein function is the Laplace exponent of a 
killed compound Poisson subordinator with Gamma distributed jumps. It 
can actually be seen as an extension to negative values for the parameter a 
in the definition of the Bernstein function considered in Example 5. It is also 
a complete Bernstein function on account of the fact that T has a completely 
monotone density. In [30] it was shown that 



A + K A + K Jo 

where p = A/(A + k). 

Thanks to the fact that is a complete Bernstein function, its conjugate 
0* may also be used to construct a scale function. From the above description 
of W one establishes in a straightforward way that k* = 0, d* = + A) 
and 

T*(a;,oo) = -f^e-^^x'-'E^APl^x'') 

A + K 
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and so 



It follows that the respective conjugate parent processes have Laplace 
exponents 

jp(e) ^K9 + xe(i- (^^] ] and tij*{e) 



Moreover, the parent process associated with (f) drifts to oo, has no Gaussian 
component and jump measure satisfying 

n(-oo, -x) = ^^x'^-^e-^'', for x > 0. 

The conjugate parent process oscillates, has a Gaussian coefficient a* — 
■\j2/{\ + k) and its jump measure satisfies 



J X 



n*(-oo, -y)dy = J^e-^-x''-'E,,,ifry''x^), x > 0. 



Note that although it has been assumed that k > 0, since both scale 
functions are continuous in their parameters through their Laplace transform, 
it follows that one may simply take limits as k J, to include k = in the 
parameter ensemble. In the latter case, the parent process will oscillate as 
opposed to drifting to oo. 

Example 7 

For q; e (0, 1], A > 0, now we take 

= Alog(l + r), ^>0. 



In the case that a = 1, is the Laplace exponent of a Gamma subordinator. 
In the case that a e (0, 1) this is the Laplace exponent of an a-subordinator, 
subordinated by a Gamma subordinator. A subordinator with (f) as Laplace 
exponent is usually called Linnik subordinator. For this Laplace exponent 
one may show that k — d — and 

T(dx) = aX ^"'^^~^^ dx, x>0. 

X 
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Note that is a complete Bernstein function as soon as one recalls that 
Ea,i(— x) (cf. [52]) and 1/x are completely monotone and that the product of 
completely monotone functions is completely monotone. The scale function 
associated with is unknown. None the less, from the discussion in Section 
[5] we see that the conjugate 0* has a Levy measure T* which is absolutely 
continuous with non-increasing density. Hence it follows from Corollary [1] 
that there exists a spectrally negative Levy process with Laplace exponent 

whose associated scale function is 

W*{x) = ax£ dy, x>0. 

It was shown in Hubalek and Kyprianou [30] that for the case a = 1 

Wix) = xj^e~yy^ Yl^dt^dy, x>0. 

It follows from the latter that d* = 0. Since (f)*{9) = 6'/Alog(l + 9), an easy 
limit as 6^ I shows that k* = 1 and hence from the expression given for W 

T*{x,oo) = Xe''' -—dt-1, x>0. 

Jo ^(t) 

It now follows that the conjugate parent process drifts to cxo, has no Gaussian 
component and, formally, Theorem [T] tells us that 11* satisfies 

oo /-oo ^i— 1 



U*{-oo,-y)dy = Xe I Y{t)^^' ^ ^ °' 



Example 8 

This example is built again on discussion found in [58]. Suppose we take 

m = log{(l + 9) + y/(l + 9y-l}, 9 > 0. 
This is a Bernstein function which has k = d = and 

T{dx) = Io{x)dx, X > 0, 
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where is the modified Bessel function of index u. It is known that the 
density of T is the Laplace transform of 



^{2z - z^)-^/^l^,^(^^^2)}dz^ dy, y>0, 



and hence is a complete Bernstein function. Following the reasoning in the 
previous example, we have the existence of an oscillating spectrally negative 
Levy process with Laplace exponent 



log{(l + ^) + v/(l + ^)2-l} 
for > 0, whose scale function is given by 

W*{x) = ^ 1^ ^i^(^z)dz'^ dy, x>0. 

In [HBj one also finds that the transition density of the subordinator, H, 
associated with is known and specifically is given by 

F(Ht e dx) = -/t(x)e-^ 

X 

and hence from ^ it follows that 

p-y ( r°° ^ 
W{x) = tk{y)dt\dy. 

Clearly we may now say that d* = and moreover from the fact that (t)*{6) = 
6/(j){6), taking limits as | tells us also that k* = 0. Moreover the 
expression for W given above also gives us 

-X roo 

T*(x,oo) = / tlt{x)dt, x>0. 

X Jo 

We see then that the conjugate parent process oscillates, has no Gaussian 
component and, formally. Theorem [U tells us that 

/ U* {-oo, -y)dy = / tlt{x)dt, x > 0. 

Jx X Jo 



35 



Acknowledgments 



We are grateful to Reinming Song and Zoran Vondracek for bringing to our 
attention several bibliographical references and for insightful comments that 
allowed us to improve the presentation of this work. This work was carried 
out over various periods for which support was provided by the EPSRC 
grants EP/D045460/1 and EP/C500229/1, and the CONCyTEG grant 06- 
02-K117-35. 

References 

[1] Abate, J. and Whitt, W. (1999) Explicit M/G/l waiting-time distribu- 
tions for a class of long-tail service-time distributions. Oper. Res. Lett. 
25, 25-31. 

[2] Avram, F., Chan, T. and Usabel, M. (2002) On the valuation of constant 
barrier options under spectrally one-sided exponential Levy models and 
Carr approximation for American puts. Stock. Proc. AppL, 100, 75-07. 

[3] Avram, F. Kyprianou, A.E. and Pistorius. M.R. (2004) Exit problems 
for spectrally negative Levy processes and applications to (Canadized) 
Russian options. Ann. Appl. Probab. 14, 215-38. 

[4] Avram, F., Palmowski, Z. and Pistorius, M.R. (2007) On the optimal 
dividend problem for a spectrally negative Levy process. Ann. Appl. 
Probab., 17, 156-80. 

[5] Bekker, R., Boxma, O. and Kella, O. (2007) Queues with delays in two- 
state strategies and Levy input. EURANDOM report 2007-08. 

[6] Berg, C. and Forst, G. (1975) Potential theory on locally compact 
abelian groups. Ergebnisse der Mathematik und ihrer Grenzgebiete, 
Band 87. Springer- Verlag. 

[7] Bertoin, J. (1996) Levy processes Cambridge University Press. 

[8] Bertoin, J. (1997) Exponential decay and ergodicity of completely asym- 
metric Lvy processes in a finite interval. Ann. Appl. Probab. 7, 156-69. 



36 



[9] Bertoin, J. (1997) Regenerative embedding of Markov sets. Probab. The- 
ory Related Fields 108, no. 4, 559-571. 

[10] Bertoin, J., Roynette, B. and Yor, M. (2004) Some connections between 
(sub)-critical branching mechanisms and Bernstein functions. Preprint 
Available via http:/ /arxiv.org/abs/math/0412322. 

[11] Bingham, N.H. (1976) Continuous branching processes and spectral pos- 
itivity. Stock. Proc. Appl. 4, 217-42. 

[12] Bondesson, L. (1992). Generalized gamma convolutions and related 
classes of distributions and densities. Lecture Notes in Statistics, 76. 
Springer- Verlag, New York. 

[13] Boxma, O.J. and Cohen, J.W. (1998) The M/G/1 queue with heavy- 
tailed service-time distribution. IEEE J. Sel. Areas Commun. 16, 
749.763. 

[14] Chaumont, L. (1994) Sur certains processus de Levy conditionnes a 
rester positifs. Stock. Stock. Rep. 47, 1-24 

[15] Chaumont, L. (1996) Conditionings and path decompositions for Levy 
processes. Stock. Proc. Appl. 64, 39-54 

[16] Chaumont, L., Kyprianou, A.E. and Pardo, J.C. (2007) Some ex- 
plicit identities associated with positive self-similar Markov processes. 
Preprint. 

[17] Chiu, S.N. and Yin, C. (2005) Passage times for a spectrally negative 
Levy process with applications to risk theory. Bernoulli. 11, 511-522. 

[18] Donati-Martin, C. and Yor, M. (2005) Further examples of explicit 
Krein representations of certain subordinators. Preprint available via: 
http://arxiv.org/abs/math.PR/0509041. 

[19] Doney, R.A. (1991) Hitting probabilities for spectrally positive Levy 
processes. J. London Matk. Soc. 44, 566-576. 

[20] Doney, R.A. and Kyprianou, A.E. (2006) Overshoots and undershoots 
of Levy processes. Ann. Appl. Probab. 16, 91-106. 



37 



[21] Doney, R.A. (2005) Some excursion calculations for spectrally one-sided 
Levy processes. Seminaire de Probabilites XXXVIII, 5-15. 

[22] Doney, R.A. (2007) Fluctuation theory for Levy processes. In Ecole d'ete 
de probabilites de Saint-Flour, XXXV 2005, volume 1897 of Lecture 
Notes in Math. Springer, Berlin. 

[23] Dube, P., Guillemin, F. and Mazumdar, R.R. (2004) Scale functions of 
Levy processes and busy periods of finite-capacity M/GI/1 queues. J. 
Appl. Probab. 41, 1145-1156. 

[24] Emery, D.J. (1973) Exit problem for a spectrally positive process. Adv. 
Appl. Probab. 5, 498-520. 

[25] Furrer, H. (1998) Risk processes perturbed by a-stable Levy motion. 
Scandinavian. Actuarial Journal. 1, 59-74 

[26] Hawkes, J. (1975) On the potential theory of subordinators. Z. W. 33 
2, 113 132. 

[27] Hawkes, J. (1976) Intersections of Markov random sets. Z. W. 37 3, 
243-251. 

[28] Hirsch, F. (1975) Families d'operateurs potentiels. Ann. Inst. Fourier 
(Grenoble) 25, no. 3-4, xxii, 263-288. 

[29] Hilberink, B. and Rogers, L.C.G. (2002) Optimal capital structure and 
endogenous default. Finance and Stochastics. 6, 237-263. 

[30] Hubalek, F. and Kyprianou, A.E. (2007) Old and new examples of scale 
functions for spectrally negative Levy processes. Preprint. 

[31] Huillet, T. (2003) Energy cascades as branching processes with emphasis 
on Neveu's approach to Derrida's random energy model. Adv. in Appl. 
Probab. 35, no. 2, 477-503. 

[32] Jacob, N. Pseudo differential operators and Markov processes. Vol. I 
Fourier Analysis and semigroups. Imperial College Press. 

[33] Kingman, J. F. C. (1967) Markov transition probabilities. II. Completely 
monotonic functions. Z. W. 9, 1-9. 



38 



[34] Kliippelberg, C. and Kyprianou, A.E. (2006) On extreme ruinous be- 
haviour of Levy insurance risk processes. J. Appl. Probab. 43, 594-598. 

[35] Kliippelberg C, Kyprianou, A.E. and Mailer. R. A. (2004) Ruin prob- 
abilities and overshoots for general Levy insurance risk processes. Ann. 
Appl. Probab. 14, 1766-1801. 

[36] Korolyuk. V.S. (1974) Boundary problems for a compound Poisson pro- 
cess. Theory Probab. Appl. 19, 1-14. 

[37] Korolyuk, V.S. (1975) On ruin problems for a compound Poisson pro- 
cess. Theory Probab. Appl. 20, 374-376. 

[38] Krell, N. (2007) Multifractal spectra and precise rates of decay in ho- 
mogeneous fragmentations. Preprint. 

[39] Kyprianou, A.E. (2006) Introductory lectures on fluctuations of Levy 
processes with applications. Universitext. Springer- Verlag. 

[40] Kyprianou, A.E. and Palmowski, Z. (2006) Quasi-stationary distribu- 
tions for Levy processes. Bernoulli. 12, 571-581. 

[41] Kyprianou, A.E. and Palmowski, Z. (2007) . Distributional study of De 
Finetti's dividend problem for a general Levy insurance risk process. J. 
Appl. Probab. 44, 428-443. 

[42] Kyprianou, A.E. and Surya, B.A. (2007) Principles of smooth and con- 
tinuous fit in the determination of endogenous bankruptcy levels. Fi- 
nance and Stochastics. 11, 131-152. 

[43] Lambert, A. (2000) Completely asymmetric Levy processes confined in 
a finite interval. Ann. Inst. Henri Poincare. 36, 251-274. 

[44] Lambert, A. (2007) Quasi-stationary distributions and the continuous- 
state branching process conditioned to be never extinct. Elec. J. Probab. 
12, 420-446. 

[45] James, L.F., Roynette, B. and Yor, M. (2007) Generalized Gamma Con- 
volutions, Dirichlet means, Thorin measures, with explicit examples. 
Preprint available via: |http: / / arxiv.org/ abs/0708.3932| 



39 



[46] Loeffen, R. (2007) On optimality of the barrier strategy in De Finetti's 
dividend problem for spectrally negative Levy processes. Preprint. 

[47] Pistorius, M.R. (2003) On doubly reflected completely asymmetric Levy 
processes. Stock. Proc. Appl. 107, 131-143. 

[48] Pistorius, M.R. (2004) On exit and ergodicity of the spectrally one-sided 
Levy process reflected at its infimum. J. Theor. Probab. 17, 183-220. 

[49] Pistorius, M.R. (2005). A potential-theoretical review of some exit prob- 
lems of spectrally negative Levy processes. Seminaire de Probabilites 
XXXVIII, 30-41. 

[50] Pistorius, M.R. (2006) An excursion-theoretical approach to some 
boundary crossing problems and the Skorokhod embedding for reflected 
Levy processes. Seminaire de Probabilites XXXIX, 287-307. 

[51] Renaud, J-F and Zhou, X. (2007) Distribution of the dividend payments 
in a general Levy risk model. J. Appl. Probab. 44, 420-427. 

[52] Pollard, H. (1948) The completely monotonic character of the Mittag- 
Leffler function {-x). Bull. Amer. Math. Soc. 54, 1115-1116. 

[53] Rogers, L.C.G. (1990) The two-sided exit problem for spectrally positive 
Levy processes. Adv. Appl. Probab. 22, 486-487. 

[54] Rosihski, J. (2007) Tempering stable processes. Stock. Proc. Appl. 117, 
6, 677-707. 

[55] Sato, K.I. (1999) Levy processes and infinitely divisible distributions. 
Cambridge University Press. 

[56] Schneider, W. R. (1996) Completely monotone generalized Mittag- 
Leffler functions. Exposition. Math. 14, no. 1, 3-16. 

[57] Song, R.; Vondracek, Z. (2006) Potential theory of special subordinators 
and subordinate killed stable processes. J. Theoret. Probab. 19, no. 4, 
817-847. 

[58] Song, R. and Vondracek, Z. (2007) Potential theory of 
subordinate Brownian motion. Lecture notes available from 
http : //www . math . uiuc . edu/ rsong/ptsbm . pdf . 



40 



[59] Song, R. and Vondracek, Z. (2007) Some remarks on special subordina- 
tors. To appear in Rocky Mountain J. of Math. 

[60] Suprun, V.N. (1976) The ruin problem and the resolvent of a killed 
independent increments process. Ukrainian Math. J. 28, 53-61. 

[61] Takacs, L. (1966) Combinatorial methods in the theory of stochastic pro- 
cesses. John Wiley &: Sons Inc., New York. 

[62] Thorin, O. (1977). On the infinite divisibility of the lognormal distribu- 
tion. Scand. Acinar. J. 3, 121-148. 

[63] Vigon, V. (2002) Votre Levy rampe-t-il? J. London. Math. Soc. 65, 
243-256. 

[64] Zolotarev, V.M. (1964) The moment of first passage of a level and the 
behaviour at infinity of a class of processes with independent increments. 
Teor. Verojatnost. i Primenen. 9, 724-733. 

Department of Mathematical Sciences 

The University of Bath 

Claverton Down 

Bath BA2 7AY 

UK. 

email: a.kyprianou@bath.ac.uk 

Centro de Investigacion en Matematicas A.C. 

Calle Jalisco s/n 

Col. Valenciana 

CP. 36240 Guanajuato, Gto. 

Mexico 

email: rivero@cimat.mx 



41 



